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| E ABSTRACT 


A theoretical study of the use of a direct reception Nd:YAG laser radar 





for target detection is reported. This work builds from a mathematical system 
model which incorporates the statistical effects of propagation through atmos- 
pheric turbulence, target speckle and glint, and receiver noise; both photo- 

; multiplier and avalanche photodiode receivers are considered. Results are 
presented for the receiver signal-to-noise ratio which show, explicitly, the 
deleterious effects of atmospheric scintillation on system performance. The 


structure and performance of the optimum likelihood-ratio processor for Single- 


pulse target detection are analyzed. 
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I. Introduction 

In a previous report [1] we examined the performance of a heterodyne- 
reception CO, laser radar in both imaging and target detection scenarios. 

The analysis contained therein was carried out only for clear-weather propa- 
gation conditions, but did incorporate realistic statistical models for 
propagation through atmospheric turbulence, target speckle and glint, and 
heterodyne-reception shot noise. 

As discussed in [2], direct-reception Nd:YAG laser radars are candidates 
for pulsed-laser ranging and projectile-tracking systems. The purpose of 
the present report is to develop a body of target detection results for a 
Nd:YAG laser system to facilitate quantitative performance evaluation in gen- 
eral and comparison with a CO, laser system in particular. In Section II, 
we shall present a mathematical system model for the Nd:YAG radar. This 
model, which draws heavily upon the results of [1], will include the effects 
of turbulence, target interaction, and receiver noise; both photomultiplier 
(PMT) and avalanche photodiode (APD) receivers will be considered. In 
Section III, we shall use this system model to derive the single-pulse 
receiver signal-to-noise ratio (SNR) for PMI and APD systems. This material, 
although not directly indicative of target-detection performance, will allow 
us to quantitatively examine the effects of target speckle and atmospheric 
scintillation. Moreover, the SNR analysis provides simple quantitative com- 


parisons between PMI and APD reception Nd:YAG radars, and may be used to 


compare Nd:YAG direct reception and co, heterodyne reception radars. The 





single-pulse target detection problem will be considered in Section IV. 

It turns out that the PMI and APD receivers both yield simple likelihood- 
ratio test (LRT) processors. Unfortunately, the associated receiver operat- 
ing characteristic (ROC) results require numerical evaluation; this evaluation 
is underway, and will be the subject of a subsequent report. Finally, in 
Section V, we present a summary of the key results we have obtained. 

II. Nd:YAG Laser Radar Configuration 

For our purposes, a direct-reception Nd:YAG laser radar may be represented 
by the block diagram of Fig. 1 (cf. [1, Fig. 3]). In drawing this diagram, 
we have made the following assumptions. 

The transmitter and receiver will be taken to be co-located. The trans- 
mitter laser will be assumed to produce a periodic train of tp sec duration 
rectangular-envelope pulses with peak power Poe These pulses will be trans- 
mitted from an unobscured circular exit pupil of diameter dy, in the range 
1mm < d, <5 m. 

The target of interest will be located a distance L from the radar along 


a line-of-sight propagation path through the clear turbulent atmosphere; L- 


values of interest will lie in the range 1 km < L < 10 kn. ‘ 
The receiver entrance optics will be assumed to be an unobscured cir- ? 
cular pupil of diameter dp in the range 10 cn < dp < 1m. The photodetector 7 | 


will be either a PMT or an APD. 


Making extensive use of the development in [1], the foregoing assumptions 


yield the following mathematical system model. 
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(1) The transmitted wave has a monochromatic linearly polarized 


electric field with complex envelope 
E,(6) = (2Pp/ce,)/? £006), (11 


for p = (x, y) a two dimensional vector in the transmitter's 
exit pupil and Ep (o) a normalized (unity square integral over 
the transmitter pupil) spatial mode. 

(2) The complex envelope of the electric field arriving from the 
transmitter in the z = L plane, E, (e') for 9’ = (x', y'), sat- 
isfies 


E,(OV% (2Pp/ce,)/? £6") 


at) 


x exp (x(0',0) + jo(o’, 09), (11.2) 


where 


£406") = fab &p(6) GAL)? expljKLG +b - 87/2] (AT. 3) 


is the normalized free-space target plane field mode generated by 


transmission of Ere) , and x, $ are turbulence-induced log- 


amplitude and phase fluctuations. Equation (II.2) assumes that 


dy < Po where 
e. ona 2ys 23/5 
Py * (1.09k C. L) 
is the turbulence field coherence length (see Fig. 2). 


(3) The complex envelope of the electric field in the z = L plane 


that is directed back towards the transmitter satisfies 





(11.4) 


10” 


Po (m) 


10° 





L (m) 


Fig.2. Turbulence coherence length p,, ¥S propagation 
path length L for weak (c2 = 5 x 10-16 m-2/5), moderate 
(c = 10°14 m7 2/3), and strong (ce = 5x 10°*° m7 2/3) 


turbulence; 1.06 um wavelength has been assumed 
throughout. 





—. 


i. E,.(e") . E, (09 T (0"), (11.5) 
where 


1(6') = 1,(6') e)° + T.6") (11.6) 


is the complex-field reflection coefficient of the target. In 
Eq. (II.6), qT, is the glint reflection coefficient, ® is a 
random phase angle, and T, is the speckle reflection coefficient. 
The characterization of these quantities given in [1, Sect. III] 
will be assumed herein. 

(4) The complex envelope of the electric field (due to the target) 


in the receiver entrance pupil, Ep(o), satisfies 
Ba(6) = _f'do" E,(6") hy (@", 8), (51.7) 
where 
a ee el : =~ . =poyage 
hy (o', ©) = (jAL) ~~ exp[jkL(1 +] 9" - o|°/2L")] 
x exp [x(o', 6) + jo(o', 0)] (11.8) 
is the stochastic atmospheric Green's function. 
The statistics of the photocurrent depend on whether PMT or APD reception 
is employed. The relevant noise models for these receivers are developed 


below. 


PMT Receiver 


An idealized PMI optical receiver is shown in Fig. 3. It consists of i 
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a circular objective lens* (focal length 2, diameter dp) to collect and 
focus signal light, an optical-frequency filter (bandwidth AA in wave- 
length units) to discriminate against background light, and a photomultiplier 
(circular active region of diameter dp) to produce an electrical signal 
proportional to the incident optical power. 

For typical PMT parameters, we may assume that the photocurrent i(t) 
is a conditional Poisson point process [3] - [5]. In particular, assuming 
the receiver is illuminated by Eg (0) during the t_-sec pulse duration, plus 


P 
a broadband spatially-distributed bates field then 


r= factio, (11.9) 


0 
which is the number of photons counted during the pulse duration, is a 


Poisson random variable conditioned on knowledge of T, x, ». This Poisson 


random variable has mean value (over the shot-noise ensemble) 
<N> = favo, (11.10) 


where u(t) is the effective average photon arrival rate at time t. When 


d,, is large enough to encompass the entire angular spectrum of ER (0) we 


D 
find that [6] 


*For values of interest, reflecting rather than refracting optics must 
be oyed. The two, however, are interchangeable from a detection theory 


viewpoint. 





rf 


u(t) % 
(n/hyg) [(cey/2) fad circ (2[8|/dq) EQ)? + PJ (HE-14) 


In Eqs. (11.9), (11.11), q is the charge released per photoabsorption, n is 


the detector's quantum efficiency, and 
eiye 2 
Pi =} NAA (mdpd,/4A2) (11.12) 


gives the average background light power incident on the detector in terms 
of the background spectral radiance [7], [8]. (Note that we have neglected 
dark current, which would add another term to u(t).) 
APD Receiver 

An idealized APD optical receiver is shown in Fig. 4. Basically, it has 
the same structure as the PMI receiver of Fig. 3. There are significant 
practical differences (a solid state APD is more rugged than a vacuum tube 
PMI, etc.) as well as marked analytical differences (basically Gaussian APD 
photocurrent statistics vs. basically Poisson PMT photocurrent statistics). 

There are available a number of detailed statistical analyses of APD 
optical receivers [9] - [11]. We shall use the simple Gaussian approximation 
theory [12], [13], which is known [11] to give fairly accurate results in 
error probability calculations, and uses exact results for the signal-to- 
noise ratio. In this approach the photocurrent i(t) is taken to be a 
conditional Gaussian process. Specifically, if the receiver is illuminated 
by ER (6) during the t,-sec pulse interval, plus a broadband spatially- 
distributed background field, then conditioned on knowledge of T, x, 9, 
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the current i(t) is a Gaussian random process with mean value (over the 


receiver noise ensemble) 


<i(t)>, = eGu(t) (II.13) 
and covariance function 


<(i(t + 1) - <i(t + 1)>,)(G(t) - <i(t)>,)>, 


= [e7G2F u(t) + 2kT/R] 6(1). (11.14) 


in Eqs. (11.13), (11.14), u(t) is the effective average photon arrival rate 
given by (II.11), e is electron charge, G is the average APD current gain, 


1/2 is the excess noise factor that is due to the randomness of the 


FXG 
APD gain mechanism, and 2kT/R is the Johnson (thermal) noise spectrum 
associated with the load resistance R. (Note that in this model we are 
neglecting the noise figure of the electrical preamplifier that follows 
the photodetector.) 
III. Signal-To-Noise Ratio Analysis 

Consider the use of the radar configuration described in Section II 


to observe a target at range L. As in [1, Sect. IV], let us assume that 


a single pulse is transmitted with spatial mode 


Eq (6) = (4/ndy?)'/? exp (j2nFp + 5), for |8| <dy/2 (IIT. 


i.e., the transmitter radiates a collimated beam which propagates in the 
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direction of the unit vector 
: = 2 1/2 : 
ip = (2mE, + (k? - |2nE|2)/? Lak 


where i, is a +z directed unit vector. For dy? /AL <1, dp < Po» the target- 


plane field E, (0') therefore satisfies (11.2) with 


&, 06") & (nd_2/ay!/? Gat)? expGikL + jk [6"|? /2L) 
x Jy (mdp|o' - ALE,|/AL)/(ndp|o" - ALE,|/2AL) . (111.2) 


Thus, the radar illuminates a nominal AL/d> diameter region of the target 
convered on the point p' = ALE. 

For realistic projectile-tracking scenarios, the target of interest will 
lie entirely within this AL/dp diameter region. Hence, in computing Ep» the 
target-return field, from (II.5) - (11.8) and (III.2) we may use 
&, (ALE,) exp (j2nf, + (o' - ALE,)) in place of &,(p'). Furthermore, it will 
typically be true that d)/L (the angle subtended in the z = L plane by the 
receiver entrance optics) will be less than the angular beamwidth of the tar- 
get return. This enables us to use p = 0 in the Fresnel diffraction kernel 
in (II.8) when we calculate ER: Finally, assuming (as in [1, Sect. IV]) 


that the target reflection coefficient may be treated in the single-glint 


approximation, we find that 








Ep(6) % (2Pq/ce,)!/? &, ALE,) GAL)? exp(jkL - j4maL |Ep|7) 


x fab" (1,6) explx@", © + x, 0) + 5G", 0) + 66%, 6))] 
+ T,(6") explx(0",, 0) + x(6",, 0) + J(0(6',, 0) * o(0",, 6) + 8)1) 
x exp (j4nf, eo"), (1II.3) 


is the target-return field in the z = 0 plane. Using (III.3) with the 
statistics given in [1] for T., xX, %, 9 will enable us to evaluate the re- 
ceiver signal-to-noise ratio for PMT and APD systems, as shown below. 
III.1 SNR for Photomultiplier Reception 

It may be seen from (II.11) that the conditional rate function u(t) 
is time-independent, for all target and turbulence conditions, during the 
pulse interval 0 < t < ty As a result the interval count N given by 
Eq. (II.9) is a sufficient statistic for estimation of target parameters 
based on observation of {i(t) :O0<t< to! {14]. Moreover, with the 


definitions 
n, = (cegnt,/2hv.) ff db circ (2|51/d,) [ER (6) 17, (IIT.4) 


and 


my, = nPht /hvy, (III.5) 


the analysis of [6] tells us that the signal-to-noise ratio 
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(<N>-n,)? 
SNR = (III.6) 
Var (N) 
obeys 
<n_> 
. (111.7) 


1 + Var (n,)/< a n,/< n,> 


where angular brackets denote averaging over the shot-noise, turbulence, and 
target ensembles. 

Physically Eq. (III.7) is directly analogous to the heterodyne-receiver 
SNR [1, Eq. (IV.7)]. The quantities n, and nh, are,respectively, the average 
number of signal counts (conditioned on knowledge of T, x, and $) and the 
average number of background counts during the pulse interval. The numerator 
in (III.7) is the quantum-limited (signal shot-noise limited) signal-to-noise 
ratio that prevails in the absence of target-return fluctuations and background- 
light shot noise. The denominator in (III.7) measures the degradation from 
quantum-limited SNR that is due to the presence of target-return fluctuations 
and/or background-light shot noise. To emphasize this interpretation we shall 
rewrite (III.7) as ; ; 

SNR 


SNR © ——____4_ (IIT.8) eT 


i SNR,/SNRe ap + (SBR)” 


where SNR, = <n,> is the quantum-limited signal-to-noise ratio, 





oer ce 


SNRe ap = <n, >*/var(n,) (III.9) 


is a saturation signal-to-noise ratio, and SBR = <n, >/m, is the signal-to- 


background ratio. When SBR > 10, (III.8) reduces to the universal curve 


SNR_/S 

aR age ee . (111.10) 
SNRgyp 1 + SNR|/SNRoqr 

as shown in Fig. 5, increasing SNR, beyond SNRo qt in this regime does not 

appreciably increase signal-to-noise ratio. 

To complete the SNR analysis of the PMI receiver, we need only find 
<n, > and Var (n,) so that SNRo ap and SBR may be evaluated. It is easy to 
calculate, from (III.3), (III.4) and [1], that 


<8? = 4h, s + <n. >, (III.11) 


gives the unconditioned average signal count <n _ in terms of its glint 


contribution 


<ng>, = (nP yt /hvgL”) lé, ALE,) |? Ap (111.12) 


xe" (Ay Eas © 1) f a cire (2|p|/d,) exp (4C,, (0, e)), 
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and its speckle contribution 
<n. >. = (nP,t_/hv L?) lé (ALE) |? Ay (111.13) 
ss Tp ‘o =t T ci 


x TALE,) f dp circ (2|p|/dp) exp (4C, (0, 6)) 


In Eqs. (111.12), (111.13), Ar is the target's area, p' ‘is the target's 

° bidirectional reflectance for a specular reflection of the illuminator beam 
back towards the radar, and i, is the average intensity diffuse-reflection 
coefficient of the target. Both of these equations are merely notational 
translations of the familiar monostatic radar equation. With somewhat more 


effort it can be shown that for apure-glint target 


Z 
a 2 Gg, | 
Var(n.), = <n, > eile x - 1), (111.14) 
and for a speckle target 
ze 160° oO oO 0 | 
Var(n.) . = <n [c"" (ex - 1) (0° +1)/0" +1/0°] (111.15) 
where oy is the log-amplitude variance (Fig. 6), 


ot = Uf db, fab, cire (2|8,|/4g) circ (2|5|/4p) 
x exp (AC, (0, 54) + 4C,(6, 82) 


x (exp(4oy + AC, (0, 64) + AC, (0, 82) + AC (G, 61 - 69)) - 1] 


te 2 
[Of db circ (2/81 /dg) exp (4C, (0, 0)))7 (e¥%x - 19] (111.16) 
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Fig.6. Log-amplitude variance of vs propagation path length 
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L for weak (Cc, = 5x 10 m ), moderate (Cc, = 10 m 

-13 m 2/3 
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) turbulence; 1.06 um wave- 


and strong (c? = 5x 10 
length is assumed and the weak perturbation theory is employed. 





and 
ot = Of db, f db, cire (215, |/dg) circ (2lB_I/dp) fa, { 5, 
x exp(4C, (0, 64) + AC, 52)) > 
x[exp(4C, (Oy! = By", D) + AC Cy! - By", By) + AC, (By ~ B's ~ 8p) 
+ AC, (By! = B's By -Bg)) - 1) 


2 
[Ag f di cire (215 1/dp) exp(4c, (0, 6)))? (el - 19) (111.17) 


are aperture-averaging factors, 


and 


p41 + nd? A,/4(AL)’, (III.18) 


is the number of free-space degrees of freedom in the diffuse-target observa- 
tion*. A result can also be obtained for the general form of Var(n,), but it 
will not be presented here; when there is no turbulence, however, it is easily 
demonstrated that Var (n,) = Var(n,),. We shall consider specific numerical 
examples of the PMI receiver signal-to-noise ratio after we develop the 
corresponding theory for the APD receiver. 
III.2 SNR for Avalanche Photodiode Reception 

Because i(t), conditioned on knowledge of T, x, $, consists of a time- 


independent mean for 0 < t < % embedded in an additive white Gaussian noise, 


*Equation (III.17)is derived under some approximations that are discussed in 
Section IV. Moreover, some care must be exercised in using (III.16), (III.17), 
as ote ae of scintillation will be encountered in the scenarios of interest 
(see Fig. 6). 
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the matched-filter output 


N = fa (eG) ti¢t) (111.19) 


is a sufficient statistic for estimation of target parameters based on observ- 
ation of {i(t) :0<t< to {15}. Defining, therefore, a signal-to-noise 
ratio for APD reception by 
(<N'> - n,)? 
SNR = ———————_- , (III. 20) 
Var (N') 


(11.13), (11.14) may be used to show that 


< > 
The 


SNR = pHTLs21) 


Cet a a Re Ee TS EPRI ROO UG ee 
F(L + n,/<n,>) + Var(n,)/<n, > + 2kTt /Re'G'<n, > 


where My, <n,>, Var (n,) are as in Section III.1. 

Equation (I11I.21) admits to an interpretation similar to that of (III.7). 
In particular, the numerator in (111.21) is the quantum-limited SNR that pre- 
vails when target-return fluctuations, background-light shot noise, APD gain- 
fluctuation excess noise, and load resistor thermal noise may be neglected 
compared to signal-light shot noise; the denominator in (III.21) quantifies 
the SNR degradations that result from the foregoing phenomena. To emphasize 


this interpretation, we rewrite (III.21) as follows 


SNR ‘ 


SNR = ———______1___,——___,, , (III.22) 
+ (STR) 


'* SNR, /SNReat + F(SBR)~ 








where SNR,» SNR ar and SBR are defined as in Section III.1, and 


STR = <n, > e"G°R/2kTt,, is the signal-to-thermal noise ratio. When 


SBR > 10, STR > 10/F, (III.22) reduces to 


SNR, /S 
SNR ~ a ae (III. 23) 
AT q AT 


; %, . ‘ ; 2 
if follows that SNR ¥ SNR¢ ap in this regime when SNR, ae SNRo ap (see Fig. 5). 
The relative merits of PMI and APD reception, as measured by a signal- 


to-noise comparison, can be qualitatively discerned from (III.8), (III.22) 


and 
SNR. APD = SNR a par Napp/ “pwr? (III.24) 
(SNRo ap) app = (SNRg ap) pyre (III. 25) 
SBRapp = SBRpwr- (III. 26) 
In particular, at 1.06 um wavelength por values range from 10> to 10-2, 


whereas "app values range from 0.1 to 0.5. Thus, the APD enjoys a substantial 
advantage over the PMI in SNR, Part of this advantage is immediately lost 
to gain fluctuations, viz. SNRapp < SNR,/F follows directly from (II1I.22), 
and F ¥ 10 for G = 100. Moreover, at very high SNR, values both the PMT and 
APD receivers will have signal-to-noise ratios limited, by target-return 
fluctuations, to the same saturation value. Finally, at low SNR, values the 


PMI receiver may vastly outperform the APD receiver, because of the thermal 


noise that is present only in the latter system. Detailed examples appear in 


the next section. 





t III.3 SNR Examples 
To illustrate the behavior of (III.8) and (III.22), we shall consider 
a series of examples analogous to those of [1, Sect. IV]. 
Case 1 No Turbulence, Specular Target 


In the absence of turbulence and target speckle, we find that, 


SNR SAT =o, aa (TII.27) : 

and : 

SNR? 4 = Ppt, (tdydp/4AL) “Ay | 

xp'(A3 Eps -E,)/hvgb.’, (111.28) | 

where the superscript "o'' denotes free-space propagation, and the subscript 
"g" denotes glint target. In Fig. 7 we have plotted the quantum-limited 

signal-to-noise ratio (111.28) vs. path length for PMT and APD reception, 


using the parameter values given in Table 1. As in [1], we include in Fig. 7 
SNR, curves given by (III.28) times exp(-2oL), for a representative value 
of the extinction coefficient a. Using the parameters given in Table 1, we 


find that 


0.206, for PMT reception 
"p 10.3, 





for APD reception 


, : and zKTt,/Re"G” = 3.23 x 10°. Combining these results with Fig. 7 gives the 


: 


SNR curves shown in Fig. 8. We see that the PMI receiver, in this case, is 
basically unaffected by background-light shot noise except when a = 0.41 kn! 
and L > 8 km, whereas the APD receiver performance shows the effect of gain- 


22 








* SMa ans an I wicca OF is ; a <7 ~ 
+ NE ae so nae Sees 5 i tak ede a IY (Ona MS Beit Par ae 
= VE a ih Ps ee ay eR > anc ee 


» 





40 


30 


SNR Gq (dB) 


20 N 
\ PMT, a=0.0 


\ APD, a=0.41 km"! 


-20 PMT, 220.41 km"! 


T 2 3 4 5 6 7 8 9 10 
L (km) 


Fig.7. Free-space propagation quantum limited signal-to- 
noise ratio vs path length L for a glint target; Eq.(III.28) 
is used with the parameter values from Table l. 
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Fig.8. Free-space propagation glint target signal-to-noise 


ratio vs path length L; Eqs.(III.8), (111.22), (111.27), and 
(III.28) are used with the parameter values from Table 1. 











TABLE I 


















PARAMETER VALUES FOR SNR’ oq EVALUATION 


Transmitter Parameters 


peak power 
wavelength 
photon energy 
pupil diameter 
pulse duration 


Channel Parameters 


path length 
extinction coefficients (see [2]) 


background spectral radiance 


Reciever Parameters 


pupil diameter 

detector quantum efficiencies 
load resistance 

mean APD gain 

load temperature 

optical filter bandwidth 
receiver field-of-view cone angle 


Target Parameters 
area 
bidirectional reflectance 


ae 
P, =.10° W 
A = 1.06 um 
a -19 
hy, = 1.87 x 107°” J 
d, = 5 mm 
t. = 50 nsec 
P 
L 
a = 0.0 
a = 0.41 km! (1.8 dB/kn) 
N, = 10 W/m? SR um 


G = 100 

T = 300° K 
‘3d 

AX = 10 * um 


dy/ 2 = 1074 


Ay * 10°? m2 


ot = 107) 


oy 


Eo ss 


fluctuation excess noise, and is strongly affected by thermal noise when 
a = 0.41 km? and L > 5 kn. 
Case 2 No Turbulence, Diffuse Target 


For free-space propagation and a speckle target we find that 
SNR® cap = 0°, (III.29) 


and 
Oise 2 
SNR eq nPrt, (mdpdp/4AL) Ar 


x T, (ALE) /hv,L”, (IIT. 30) 


where the subscript ''s" denotes speckle target. The Ny» 2kTt /Re°G? values 
for Case 1 apply here as well, thus, assuming a. = 0.1 (arbitrarily) we get 
SNR? 0 = SNR’ og and the signal-to-noise ratio curves of Fig. 9. Comparing 
Figs. 8 and 9 shows that target fluctuations are the dominant performance 
limiting effect when a =0.0 with 1 < L < 10 kn, and when a = 0.41 kn} 
with 1 < L < 3 kn. 
Case 3 Turbulence, Specular Target 

For a specular target viewed through atmospheric turbulence, we can 


show that 


2 
SNR car /(e*y - 1), (IIT.31) 
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Fig.9. Free-space propagation speckle target signal-to-noise 


i ratio vs path length L; Eqs.(III.8), (III.22), (III.29) and 
: (111.30) are used with the parameters from Table 1 and ‘ = 0.1. 





: and SNR g x SNR og? when oy < 0.5 and dp?/AL >> 1*. Using these results 
and the parameter values from Table 1, we get the SNR curves of Fig. 10. In 

| Fig. 10, we have assumed a modest turbulence strength «7 =5x10)5 no! sy, 
and have cut off the plots at the onset of saturated scintillation. Comparison 
of Figs. 8 and 10 reveals the severe SNR limit set by scintillation on a glint 
target. 
Case 4 Turbulence, Diffuse Target 


For a speckle target viewed through atmospheric turbulence, we can show 


that 
o” III.32) 
SNR car * 0 doe eo 
1+ (0 + 1)(e™yx - 1)/01 + AL) 
and SNR.o ~ SNR® og? when 0 < 0.5 and dp2/AL >>1%**, Assuming T, = 0.1, the 


*Physically, the target-return field is subject to both transmitter-to-target 
path scintillation and target-to-receiver path scintillation, and ¢' is the 
relevant aperture-averaging factor for these fluctuation terms when the tar- 
get is specular. Because dy < Qo» there is no aperture averaging of the 
transmitter-to-target scintillation in ¢', but, because dp>> dy, there is, 
in general, aperture-averaging of the target-to-receiver scintillation in ¢'. 
Indeed, when of < 0.5 so that the weak perturbation theory is valid, 
dp2/AL >>1 impties that target-to-receiver path scintillation is completely 
averaged out. 


*kFor a speckle target, the appropriate aperture averaging factor is ¢". 
Physically, c'' includes the aperture averaging of both the transmitter-to- 
target and the target-to-receiver fluctuations that occurs in the target plane, 
in addition to the aperture averaging of the target-to-receiver fluctuations 
that occurs in the recelyer pupil. When o¢% < 0.5, so that the weak perturba- 
tion theory is valid, 4 /xL >> 1 implies Xthat the target-to-receiver path 
scintillation is completely averaged out. 
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Fig.10. Atmospheric propagation glint target signal-to-noise 
ratio vs path length L; Eqs.(III.8), (III.22), and (III.31) 
are used with the parameter values from Table 1. 





parameters of Table 1, and = =5x 10745 m7 2/ 3 leads to the SNR curves 
of Fig. 11. Comparison of Figs. 9 and 11 shows that scintillation costs 
about 10 dB of signal-to-noise ratio on a speckle target whose free-space 
SNR exceeds 10 dB. Comparison of Figs. 10 and 11 shows the beneficial effect 
of target-plane aperture averaging on scintillation-limited signal-to-noise 
ratio. 
IV. Target-Detection Analysis 

Suppose there are two possible hypotheses (target present or target 
absent at a particular transverse coordinate in the plane z = L). If no 
target is present (hypothesis Hy), then the rate parameter y(t) for the 
(PMT or APD) receiver's photocurrent satisfies p(t) = nP, /hv, for 0<t< ty: 
If there is a target present (hypothesis H)) then the rate parameter for the 
photocurrent obeys (II.11) for 0 < t < ty . Our objective is to process the 
data {i(t) :0<t< ty} in such a way that the conditional probability we 
say there is a target present given there is one there is maximized, subject 
to the constraint that the conditional probability we say a target is present 
given there is no target there does not exceed a specified level. 

Because, under either hypothesis, y(t) is independent of time during the 


pulse interval, the interval count N and the matched filter output N' are 


ae Rapin 


sufficient statistics for optimally deciding between Hy and H, based on 
observation of i(t) during 0 < t < to for the PMT and APD receivers, respectively. 


The structure and performance of the optimum PMT and APD receivers will be 


developed in this section. 
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Atmospheric propagation speckle target signal-to-noise 
ratio; Eqs.(III.8), (III.22), and (III.32) are used with the 
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IV.1 Target Detection with a PMT Receiver 

To maximize the probability of detection (conditional probability you 
say Hy given H, is true) at a given false alarm probability (conditional 
probability you say H, given H) is true) based on knowledge of N, we must 
use a randomized likelihood-ratio test, because N is a discrete random 
variable[16]. The structure of this test iseasily developed from well known 
results for doubly stochastic Poisson processes [17], [18]. 

The likelihood-ratio test (LRT) for the PMT receiver is 


say Hy 


> 
~ 


<Q + n/n)’ exp-n> oH, qv.) 
< 


say Ho 
where N is the received photon count during the pulse interval, n . and My 
are as defined in Section III, and angular brackets denote averaging over the 
turbulence and target ensembles. The threshold value n in (IV.1) will be 
determined as described later. Because the left-hand side of (IV.1) is a 
monotonically increasing function of N, it is equivalent to the simple thresh- 


old test 


en eee : 


say Hy 


> 





En. 


<a 


fred ¥ 





where y is an integer chosen as described below. Note that the processor 
(1V.2) gives the LRT regardless of the strength of the turbulence, or whether 
the target is pure glint or pure speckle. ) 

The receiver operating characteristic (ROC) for (IV.2) is a discrete set 
of points {(PL(), Pyty)) ‘y= 0, 1, 2, « « .} where Prly), Py(y) denote the 
false-alarm and detection probabilities obtained by using (IV.2) with thresh- 
old y. As shown in Fig. 12, a discrete ROC cannot provide a test which achieves 
maximum P, at arbitrary P, values {16]. Suppose, see Fig. 12, that p is the 
maximum allowed false-alarm probability and that Yp is a non-negative integer 
such that 


Pe(y)) > P> Pelyy + (IV.3) 


Using the randomized test, "say H, when N =% + 1, say Hy when N < ip and 
when N = Yp say H, with probability § and Hy with probability (1 - 8)" gives 
false-alarm probability BPE(Y,) + (1 - g) PE Yp + 1), and detection probability 
BP (y,,) * (1 > BP %) (Yp + 1); these probabilities appear as the straight line 


in Fig. 12 as 6 varies from 0 tol. It is clear from this line that choosing 
B= (P - Pry + 1))/(Pplyp) - Ppl + 10) (IV.4) 


gives a test which has false-alarm probability exactly equal to p, and a higher 


detection probability than any LRT whose false-alarm probability is less than 


‘ p; a more refined argument [16] shows that no test has a higher detection 


probability at false-alarm probability p or less. 
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., Fig.12. Illustration of a discrete random variable receiver 
operating characteristic; the open circles show {(PR(Y), 

P (y)): 0 € ¥ € 5}; the solid line is the randomized test 

per ormance for Y, = 1; the intersection of the solid line 

with the Pr =p tine gives the maximum Py consistent with Pes p. 





The foregoing analysis, and Fig. 12 in particular, shows it suffices to 
derive the ROC {(P.(y), Py(y)) >: y = 0,1, 2, .. . } to characterize the 
performance of the LRT (IV.2); randomized test performance is then found 
by connecting adjacent ROC points by straight line segments. It turns out 
to be relatively simple to deduce (Pry), Pyly)) when the target is either 
specular or diffuse. Regardless of target and turbulence conditions, N is 
a Poisson random variable with mean n, under Hy - Thus, for any target/tur- 


bulence combination we find that 


co 


Pay) = = (nt)"* (n,)” exp(-n,). (IV.5) 


n=y 


The general expression for Pyty) is 


oo 
-1 n 
Ppt =< Ee eh ig 85) er, ay) (IV.6) 
where angular brackets denote averaging over the turbulence and target 
ensembles. More explicit results are developed below. 
Case 1 Specular Target 
For a specular target, we find, from (III.3), (III.4), (III.12), that 


2w 
n, <n, >,€ ; (IV.7) 


























where 


a, Lib sre Clap ool’ D+ Gy DI 


“ : ae (IV.8) 
dp circ (2]p|/dp) exp(4C, (0, 6)) 





Because x(6', p) is a Gaussian process, and sums of real-valued lognormal 
random variables are themselves approximately lognormal random variables 
[19], [20], we have that w is a Gaussian random variable. By direct calcula- 
tion of the variance of n, it follows that w has mean -o'* and variance o'” 
where 


2 2 
gr et eee hy = 1) (Iv.9) 


and ¢' is given by (III.16). Thus, the glint-target PMT receiver detection 


probability satisfies 


2W 
Po, -f av PyW) Py” (ys <n, mes Mh) (IV.10) 


with p,(W) = (2102) "1/2 exp[-(W + o!*)2/20'7], and 


Py'(ys a, b), = E nt)" (a+ b)” exp[-(a+b)]. — (IV.11) 
n=y 


The quantity Py (ys a, b) g given by (IV.11) is the glint-target free-space 
detection probability achieved by the test (IV.2) when the average target 
return count is a and the average background count is b. Equation (IV.10) 


shows, therefore, that the presence of turbulence introduces a lognormal 


fading into the detection probability calculation. 








Case 2 Diffuse Target 

For a diffuse target, (III.3) shows that Ep (0) is a circulo-complex 
Gaussian random process conditioned on knowledge of x, ». By means of the 
equal-eigenvalue approximation, we then have the conditional counting dis- 


tribution [21]- [23]. 
PrIN=n |x, 6, H,] % 


[(A,/0)"/ (1 + n.yo)"* 9) exp[-n,/(1 + A,/0)] 


0-1 


xL, [-n,0/n, (1 + n,/0)}, (IV.12) 


where 
% 2 2 = 
n, = (nPt,,/hvob a & OLED) | T, LED) 
x fap circ (215|/dq) f a5" exp [2x(6', 0) + 2x(6',8)1, (1V.15) 


is the average number of counts due to the diffuse target conditioned on know- 
ledge of the turbulence parameters. In Eq. (IV.12) i “1 is the generalized 


Laguerre polynomial of index n and order 0-1, i.e., 


n /n+D-1 : 
ie we = £ ( (-x)*/il, (IV.14) 
i=0 i 


and D is the number of degrees of freedom of the observation [22]. Strictly 


speaking, D is a function of the turbulence parameters x and $; this makes 











subsequent averaging of (IV.12) over the turbulence ensemble extremely dif- 
ficult. We shall use the reasonable approximation [24], that D is given 
approximately by its free-space value 0°, which, for a target that lies 
entirely within the AL/d> - diameter illumination region, satisfies (III.18). 
Now, using once more the convergence properties of sums of lognormal 


random variables, we may express n, in the form 


Si 2v 
n, = <n, >,€ (IV.15) 


where v is a Gaussian random variable with mean =? and variance ott, and 


2 
eto" 16 


2 
-1l=e ce °%y - 1). (IV.16) 


We now have that the speckle-target PMT receiver detection probability is 
eo 
P.(y) = {dvp.) P-° (y; <n, >.e”, n,) (IV.17) 
pty Vv p S¥e My 75) » Bhs? : 
-0 


where p,(V) = (2nott2) “1/2 exp[-(V + ott?) 2/297], and 
‘ co Oo 
Pr? (ys a, bY, = 5 [(@/0°)"/( + a/0°)"*" ] exp[-b/(1 + a/0%)] 
n=y 


xi”) [-po%aca + 2/01, (1V.18) 


is the free-space speckle target PMI receiver detection probability when a is 
the average target-return count and b is the average background count. 
Calculations are available [22], [25] for the free-space PMT receiver 


ROC, for both pure glint and pure speckle targets; some examples will be given 
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after we analyze the ROC for APD reception. Calculations have yet to be 
. done for the PMT receiver ROC in the presence of turbulence. 
IV.2 Target Detection with an APD Receiver 
Well known results for Gaussian detection [15] permit us to show that 


the likelihood-ratio test based on the observation of N' is 


<[1 + n/(m, + 2kTt/Re?67R) }"/ 2 
x exp(-[N'- (n, + n,)17/2[(n, + nF + 2kTt,/Re"G"] 
+ (Nt - n,)7/2[n,F + aKrt »/Re?G"]}> 


say Hy 
= 
n (IV.19) 


< 
say Hy 
where angular brackets denote averaging over the target and turbulence ensembles, 
and the threshold ¥ is chosen to achieve the false alarm probability constraint 
with equality. 

It turns out that the left member of (IV.19) is a monotonically increasing 
function of N', thus (IV.19) can be reduced to the simple threshold test 
say Hy 
= 


Nt > (IV. 20) 


say Hp 








where y is chosen to satisfy the false alarm constraint with equality. This 
test is the optimm processor for all target/turbulence combinations. 

For the test (IV.20), it is easy to find the false-alarm probability 
Pr» because under Hy N' is a Gaussian random variable with mean nm and 
a 


variance age = nF + 2kTt,,/Re Thus for all target and turbulence conditions 


we find that 

‘i Pe = QI(y - n,)/o9], (IV.21) 
where Q(x) = (2m) “1/ 2 f dy exp(-y2/2) is the area under the tail of the nomal- 
ized Gaussian density.* Equation (IV.21) can be inverted numerically to obtain 
y as a function of P_, nh, and Do: 

The general result for the detection probability Py is obtained by noting 
that given T, x, 9, H)> the random variable N' is Gaussian with mean no+n, 
‘ Be ane 

and variance oy = (n, + n,) F + 2kTt /Re G“, whence 

Py = <QICy - (my + m,))/0,]>. (IV.22) 
In what follows we consider the pure glint and pure speckle limits of (IV.22). 


Case 1 Specular Target 


For a specular target, we have that A, = <n, hs from (IV.7). This 


tells us that the glint-target APD receiver detection probability obeys 
P pty g a 


fate, 00 QLCy- Cen, >,07" + m))/0, 001, (1v.23) 
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where P,,(W) = (2no' 2) “1/2 exp[-(W + ot 7) 2/207] as before, and 0,7) = 

2W ene 
(<n, >ge + n,) F + 2kTt,,/Re G”. 
Case 2 Diffuse Target 

By means of the equal-eigenvalue approximation, with the number of free- 
space degrees of freedom 0° used in place of the number of atmospheric degrees 
of freedom D, we can show that the speckle-target count parameter n. is given 
by 
a 2v 
no = <n > me’, (IV. 24) 
for v is as in Section IV.1, and m a random variable that is statistically 
independent of v with density function 
P.M) . [0° (omy™ - 1) cp? - 1)!] exp(-0%™) u™), . (IV.25) \ 


where u(M) is the unit step function. Thus, the speckle-target APD receiver 


detection probability is 
Pps = far) f ap,00 


x QL(y - (cng >, Me + n,))/o,M, V)] (Iv. 26) 


where 0,7(M, V) = (<n.>,MeV+ n,) F + 2kTt /ReG”. 


As yet ROC calculations are not available for the APD receiver, except 
when it is thermal-noise or background shot-noise limited [15], [26]. 
_IV.3 ROC Examples 


The numerical evaluation of the foregoing ROC results is presently 


underway, and will be the subject of a subsequent report. In this section, 
» 





{ we shall present some simple free-space performance results. 
Case 1 No Turbulence, Specular Target 

In the absence of turbulence, the PMI receiver glint-target detection 
probability is given by (IV.11) and the randomized test procedure (IV.3), 
(IV.4). We have plotted this detection probability in Fig. 13 for several 
false-alarm rates with n, = 0.0 andn, = 0.2. In this figure, we have used 
SNR? a in lieu of <n, >» to facilitate obtaining detection probabilities for 
the parameters given in Table 1 via Fig. 7. We see that, for the parameters 
given in Table 1, a 95% detection probability can be maintained with 
P, = 10° out to a S km path length when a = 0.41 kn"? 

The APD receiver free-space glint-target detection probability, obtained 
from (IV.23) with ot? = 0, is shown in Fig. 14 for several false-alam rates, 
assuming the parameters of Table 1. We see, from Figs. 7, 14 that the APD 
system can achieve Py = 0.95, PE = 107° almost out to L = 5 km when 
a = 0.41 kn. This is so because the 17 dB SNR, advantage of the APD 
receiver just about balances out its thermal-noise performance degradation 
Case 2 No Turbulence, Diffuse Target 

In Fig. 15 we have plotted the free-space speckle-target PMT receiver 
detection probability, assuming no background noise, for several degrees of 
freedom values. Comparison of Figs. 13 and 15 shows that when p° = 1 we 
need about an 8 db SNR”, increase to maintain Py = 0.95. On the other hand, 
for 0° = 102 speckle target and glint target performance are virtually 


identical. The latter result is a manifestation of the constant-intensity 
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i 
(low photon coherence) limit [21], [22], which is also valid in the presence 
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Fig.13. Free-space propagation glint target detection prob- 
ability for the PMT receiver vs quantum limited signal-to- 
noise ratio. 
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Fig.14. Free-space propagation glint target detection prob- 
ability for the APD receiver vs quantum limited signal-to- 
noise; the parameter values from Table 1 are used in calcu- 
lating the background and thermal noise levels. 





99.99 


1-5-9792 


99.9 


90 


80 


50 


(percent) 


° 
* 


0.01 





’ SNR¢, (48) 


Fig.15. _Free-space propagation speckle target detection 
probability vs quantum limited signal-to-noise ratio. 








of background noise. Thus, for nm, = 0.2, Fig. 13 also gives the speckle 
target detection probabilities for the parameters of Table 1, as 0° values 
of at least 70 prevail when 1<L<5 km. A similar low photon coherence limit 
exists for the APD receiver, so that for 1<L<S kn, Fig. 14 also gives the 
speckle target detection probabilities for the parameters of Table 1. 
V. Summary 

In this section we shall summarize the key results of our signal-to-noise 
ratio and target detection analyses. We shall presume familiarity with the 
notation summarized in Table 2, and we shall not explicitly state all the 
conditions required for the validity of the results. 
PMT-Receiver Statistical Model © 

The interval photon count N for the PMI receiver is a conditionally- 
Poisson random variable. Specifically, given the turbulence and target char- 


acteristics, N is Poisson-distributed with mean n, + n, where 
ng = (nPgty/hvgL”) |g (ALEp) |? exp (-2 aL) 
x Apo" (as Eps - Ep) ff db circ (2/|/dg) expl2x(og', 6) + 2x (6g, 8)] 
+ x? f a circ (2|5|/dg) | fe T,(o") explx(o", 5)* x(6', 6) 
+ 5(6(6"» 0) - (6's 6) + 4nEp + 5D1/7, (v.1) 


and 


my = mr7aaN, (mdgdy/4A2)"t /hv,. (V.2) 
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TABLE 2 
SUMMARY OF NOTATION 


Radar Parameters 


Pr 


ct 
sc 


Me 


= 


Yo 


”s - © 2 s 


AA 
d/e 
& 


transmitter peak power 

pulse duration 

transmitter optics diameter 
receiver optics diameter 

1.06 um wavelength 

photon energy at 1.06 um wavelength 
detector quantum efficiency 

APD load resistor 

APD current gain 

APD excess-noise factor 

thermal energy 

optical filter bandwidth 
receiver field-of-view cone angle 


normalized target-plane illumination beam 


Atmospheric Parameters 


turbulence-induced log-amplitude fluctuation 
turbulence-induced phase fluctuation 
path length 


log-amplitude variance 





ee 


TABLE 2 (Continued) 


oe 


a = extinction coefficient 
’ c', c’’ = aperture-averaging factors 
o, o? = aperture averaged log-amplitude variances 


N 


| 1 * background-light spectral radiance 
| 


Target Parameters 
Ay = target area 


T. = diffuse- reflection coefficient 


5s 
o' (3 fyi - £7) = target bidirectional reflectance for a glint return 
T; = average diffuse reflection coefficient 

p = free-space propagation speckle target degrees of freedom 


aes 


ms 


eat 
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APD-Receiver Statistical Model 


———_ OO 


The matched-filter output N' for the APD receiver is assumed to be 
Gaussian-distributed, given the turbulence and target characteristics, with 
ee 22 
mean n, + ny, and variance F(n, + n,) + 2kTt,/Re G*. 
Receiver Signal-to-Noise Ratio 


The receiver signal-to-noise ratios defined by 


(<N>- n,) 2 
I ¢ 2 (V.3) 
Var(N) 
and 2 
(<N'> - n,) 
1 (V.4) 
Var(N') 
can both be put into the standard form 
SNR = hes a (V.5) 
F' + SNR q/ NRap + FF (SBR) “+ + Spr (STR) 
by the following identifications: 
quantum-limited signal-to-noise ratio 
SNR, SRL A55 (V..6) 
excess-noise factor 
p=), for PMI reception (V.7) 


F, for APD reception; 


saturation signal-to-noise ratio 


SNRo at = <n, >*/var(n,) ; 








signal-to-background ratio 


SBR = <n, >/m, 3 (V.9) 
and signal-to-thermal ratio 
STR = <n, > Re7G?/2kTt,, (V.10) 
At very low SNR. values, (V.5) shows that performance is background 


q 
and/or thermal noise limited; in this regime the SNR increases as the square 
of SNR, - For SNR, values large enough to give STR > 10/F' and SBR > 10, 
(V.S) reduces to the universal curve (Fig. 5). 


Sec % rte igs (v.11) 
Increasing SNR, beyond F'SNRo ap does little, therefore, to improve receiver 
signal-to-noise ratio. 
The foregoing signal-to-noise ratio results have the following behavior 


in particular cases of interest: 


i) No Turbulence, Specular Target (Figs. 7, 8) 
o 2 
SNR og = mP rt (ndydp/4AL) Ay exp(-2aL) 


x p'(A; fF “ D/hv,L, (V.12) 


SNROSAT = © (v.13) 








ii) No Turbulence, Diffuse Target (Figs. 7, 9) 


_ 


SNRO, = nPpty (ndpdp/4AL)7 Ay exp(-2aL) 


x T,(ALE_)/hy,L?, (v.14) 
SNRE sar = 7° (V.15) 


iii) Turbulence, Specular Target (Figs. 7, 10) 


~% oO 
SNR og v SNR og? (V.16) 
SNR can © 1/(e4° a 1) (v.17) 

gSAT V x : , 
iv) Turbulence, Diffuse Target (Figs. 7, 11) 

oO 
SNR, x SNRS 4 (V.18) 
0° ) 

SNR cap & ——>———22 (v.19 
SSAT™ 1 + (0? + 1) (eX- 1/1 + AVAL) 
Likelihood-Ratio Tests 
P The PMT-receiver decision rule which maximizes the conditional probabil- 


ity you say a target is present given there is one present, subject to the 


constraint that the conditional probability you say a target is present given 


: there is no target there does not exceed p, is the randomized test: 








i) N>v¥y, + 1, say Hy 


a 
cs 


ii) N= y,, say H, with probability 6 and H) with probability (1 - 8) 
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iii) N< Yp? say Hp, where Yp is the integer which satisfies 


oc oo 


r (n!)72 mle ?> p> 3 (n!)72 nove 
ney, tr ney, +1 


™ 


for p = the maximum allowed false-alarm probability, and 


Mh Yp .% 


B= ([p- «f= (nt)? ne 
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The APD-receiver decision rule which is optimum under the foregoing 
detection probability criterion is the threshold test 


say Hy 
ka 
N' ¥ <s (V.20) 
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say Hy 
where y is chosen to satisfy 


‘ . 2n2,1/2,.. 
P = Qi(y-n,)/(n,F + 2kTt,/Re"G")* *] 


for p = the maximum allowed false-alarm probability, and Q(x) = (2m) / 2 


fe exp (-y2/2). : 








Receiver Operating Characteristics 


For the PMI-receiver, the randomized threshold test gives detection 













probability 
Py = BPp(yp) + (1-8) Pp Cy, +) (V.21) 


where Poly) is the detection probability for the threshold test 


say Hy 
z 
N Y, Y an integer. (V.22) 
< 
say Hy 
The behavior of Pyly) is as follows 
i) No Turbulence, Specular Target (Fig. 13) 
o a les eggpe 
PD Cy) = Ph Cys SNRQ,, m), (V.23) 
where Ph (y; a,b), = £ (nl)? (a +b)” exp[-(a + b)]. 
=y 
ii) No Turbulence, Diffuse Target (Fig. 15) 


Ph Gy) = Pp (ys SNREG» m), (v.24) 
where 


Po (ys a, bY, = £ [(a/D)"/(L+a/02* 7) expl-b/(a + a/°)] 
n=y 


x - 1 F-pp°/a(1 + a/0°)]. 


iii) Turbulence, Specular Target 


Poy) =f diip,(W) PD Cvs SNRyg €™, m) 5, (v.25) 
where a 

p(w) = (2not2)"/? expf- (w+ 0'7)7/20'7), 
wen OO hs x1 (el6oy ok 


iv) Turbulence, Diffuse Target 


Py(v) =f avpy(V) Ph (1s SWRg es my), (v.26) 
where 
py(V) = (2m)? expt-¢v + o'7)7/20"7) 
2 2 : 
with e407 - 1 = (eho, - 1), 


For the APD reciever, the threshold test (V.20) gives the following 
detection probabilities: 
i) No Turbulence, Specular Target (Fig. 14) 
Pp QECy - SNR, = m,)/o4], (v.27) 
. 2s 0 2.2 
with Oy F(SNR,, + ny) + 2kTt,/Re G". 
ii) No Turbulence, Diffuse Target 


PH 7 fen QC ~ SHREG M - %)/0, 00) (v.28) 


with 0,70 = F(SNRE, M+n) + 2kTt, /Re"G”, 








oO 
and pp (M) = [0°(0°M)” ~ 1/0 - 1)!] exp (-0°M) u MM. 


iii) Turbulence, Specular Target 


Py =f dip) Q LCy - SNR, eM - 3/0, (H)] 


2 a 2W Zane 
where o} (W) = F(QR,. e+ n,) + akTt /Re G”. 


iv) Turbulence, Diffuse Target 


co oo 


Py = fav fap, 0) Py 00 Q I(y- SNR Me™ - m)/o,M, VDI, 


=o 


2 ‘ av : 22 
where o} (M, V) F(SNR,, Me + n,) + 2kTt,,/Re G". 


(V.29) 


(V.30) 











ACKNOWLEDGMENTS 


In preparing this report, the author benefited from the encouragement 








of A. B. Gschwendtner and H. Kleiman, and useful technical discussions with 





R. C. Harney. Their support and guidance have made this study possible. 





” 


REFERENCES 


r~ 


1. J. H. Shapiro, "Imaging and Target Detection with a Heterodyne-Reception 
Optical Radar,'' Project Report TST-24, Lincoln Laboratory, M.1I.T. 
(13 October 1978). 


2. R. J. Becherer, "Pulsed Laser Ranging Techniques at 1.06 and 10.6 im," 
Project Report TT-8, Lincoln Laboratory, M.I.T. (19 March 1976), 


° DDC AD-A024557/1. 
3. E. V. Hoversten, "Optical Communication Theory," in Laser Handbook, 
’ F. T. Arecchi and E. O. Schulz-DuBois, Eds. (North-Holland, Amsterdam, 
1972). 


4. D. L. Snyder, Random Point Processes (Wiley, New York, 1975), Chap. 2. 


5. R. M. Gagliardi and S. Karp, Optical Communications (Wiley, New York, 
1976), Chap. 3. 


6. J. H. Shapiro, "Imaging and Optical Communication through Atmospheric 
Turbulence ,"' in Laser Beam Propagation through the Atmosphere, 
J. W. Strohbehn, Ed. (Springer-Verlag, Berlin, 1978). 

7. W. K. Pratt, Laser Communication Systems (Wiley, New York, 1969), Chap. 6. 


8. N. S. Kopeika and J. Bordogna, "Background Noise in Optical Communication 
Systems,"' Proc. IEEE 58, 1571-1577 (1970). 


9. S. D. Personick, ''New Results on Avalanche Multiplication Statistics 
with Applications to Optical Detection," Bell Syst. Tech. J. 50, 
167-189 (1971). 


10. S. D. Personick, "Statistics of a General Class of Avalanche Detectors 
with Applications to Optical Communication," Bell Syst. Tech. J. 50, 
3075-3095 (1971). 


11. S. D. Personick, P. Balaban, J. H. Bobsin and P. R. Kumar, "A Detailed 
Comparison of Four Approaches to the Calculation of the Sensitivity of 
. Optical Fiber System Receivers," IEEE Trans. Commun. COM-25, 541-548 
(1977). 


' 12. S. D. Personick, "Receiver Design for Digital Fiber Optic Communication 
i Systems ,"" Bell Syst. Tech. J. 52, 843-886 (1973). 


13. S. D. Personick, ''Receiver Design for Optical Fiber Systems," Proc. IEEE 
65, 1670-1678 (1977). 





14. 
15. 


16 


17. 


18. 
19. 


20. 


21. 


22. 


23. 


24. 


25. 


26. 


D. L. Snyder, Random Point Processes (Wiley, New York, 1975) Chap. 5. 


H. L. Van Trees, Detection, Estimation and Modulation Theo Part I 
- (Wiley, New York, 1568), Chap. a. 
. H. L. Van Trees, Detection, Estimation and Modulation Theo Part I 
’ (Wiley, New York, 1568), Chap. Z. 


R. M. Gagliardi and S. Karp, Optical Communications (Wiley, New York, 
1976), Chap. 7. 


D. L. Snyder, Random Point Processes (Wiley, New York, 1975), Chap. 6. 


R. L. Mitchell, "Permanence of the Lognormal Distribution," J. Opt. 
Soc. Am. 58, 1267-1272 (1968). 


B. K. Levitt, "Detector Statistics for Optical Communication through 
the Turbulent Atmosphere," Quart. Prog. Rept. 99, Research Lab. Electron., 
M.I.T., 114-123 (October 1970). 


R. M. Gagliardi and S. Karp, Optical Commmications (Wiley, New York, 
1976), Chaps. 2, 3. 


J. W. Goodman, ''Some Effects of Target-Induced Scintillation on Optical 
Radar Performance," Proc. IEEE 53, 1688-1700 (1965). 


L. Mandel, "Fluctuations of Photon Beams: the Distribution of Photo- 
electrons," Proc. Phys. Soc. 74, 233-242 (1959). 


J. H. Shapiro, "Normal-mode Approach to Wave Propagation in the Turbu- 
lent Atmosphere,'"' Appl. Opt. 13, 2614-2619 (1974). 


C. W. Helstrom, tum Detection and Estimation Theory (Academic Press, 
New York, 1976), pp. -217, 


D. L. Fried and R. A. Schmeltzer, "The Effect of Atmospheric Scintillation 
on an Optical Data Channel-Laser Radar and Binary Communivations," Appl. 
Opt. 6, 1729-1737 (October 1967). ° 












UNCLASSIFIED 
SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


. IONS 
REPORT DOCUMENTATION PAGE 


1, REPORT NUMBER 2. GOVT ACCESSION NO. | 3. RECIPIENT'S CATALOG NUMBER 
ESD-TR-78-290 


4 TITLE (and Subtitle) (PE OF REPORT & PERIOD COVERED 


palatine Hi}. t t. | 
[ees Detection with a Direct-Reception Whe ae Fev fo 
































Optical Radarg 


\ 


6. PERFORMING ORG. REPORT NUMBER 
Project Report TST-27 


8. SZQONTRACT OR GRANT NUMBER/(s) 


F1962 8-78-C-Ppp2 


10. PROGRAM ELEMENT, PROJECT, TASK 
* AREA & WORK UNIT NUMBERS 


Program Element No. 65705F 
















eth 


wAUTWOR(S) 


t 


f pastivey H. Shaptxo / 


. PERFORMING ORGANIZATION NAME AND ADDRES 
Lincoln Laboratory, M.I.T. 
P.O. Box 73 

Lexington, MA 02173 




































. CONTROLLING OFFICE NAME AND ADDRESS 
Air Force Systems Command, USAF 
Andrews AFB 

Washington, DC 20331 







. MONITORING AGENCY NAME & ADDRESS (if different from Controlling Office) 










Electronic Systems Division 
Hanscom AFB 
Bedford, MA 01731 







« DISTRIBUTION STATEMENT (of this Report) 


Approved for public release; distribution unlimited. 





. SUPPLEMENTARY NOTES 








None 


. KEY WORDS (Continue on reverse side if necessary and identify by block number) 

















optical radar target speckle signal-to-noise ratio 
direct reception glint optimum likelihood-ratio processor 






atmospheric turbulence atmospheric scintillation single-pulse target detection 











|. ABSTRACT (Continue on reverse side if necessary and identify by block number) - 





» A theoretical study of the use of a direct reception Nd:YAG laser radar for target detection 
is reported. This work builds from a mathematical system model which incorporates the statistical 
effects of propagation through atmospheric turbulence, target speckle and glint, and receiver noise; 
both photomultiplier and avalanche photodiode receivers are considered. Results are presented . 
for the receiver signal-to-noise ratio which show, explicitly, the deleterious effects of atmospheric 
scintillation on system performance, The structure and performance of the optimum likelihood-ratio 
processor for single-pulse target detection are analyzed, 















1yan7 1473 EDITION OF 1 NOV 65 IS OBSOLETE 





